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Nonlinear magneto-optical rotation is a sensitive technique for measuring magnetic fields. Here, 
the shot-noise- limited magnetometric sensitivity is analyzed for the case of optically-thick media and 
high light power, which has been the subject of recent experimental and theoretical investigations. 
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I. INTRODUCTION 

Resonant nonlinear magneto-optical rotation (NMOR) 
jj], H has been the subject of extensive theoretical and 
experimental studies because it provides a very sensitive 
way to measure magnetic fields (see, e.g., Ref. While 
NMOR experiments are usually carried out with vapor 
samples of moderate optical thickness (no more than ~ 2 
absorption lengths), the intense research of recent years 
on electromagnetically-induced transparency (EIT; see, 
e.g., Ref. H for a reviewVhas motivated investigation of 
optically thick media J|, g] . 

Recently, NMOR in the vicinity of the Dl and D2 
lines by an o ptic ally thick vapor of rubidium was stud- 
ied § |, |, [10], The authors of Ref. @ used a 
5 cm-long vapor cell containing 87 Rb and a laser beam 
tuned to the maximum of the NMOR spectrum near the 
Dl line (laser power was 2.5 mW, beam diameter ~ 2 
mm). They measured maximum (with respect to laser 
frequency and magnetic field) polarization rotation ip m ax 
as a function of atomic density. It was found that ip m ax 
increases essentially linearly up to n ss 3.5 x 10 12 cm~ 3 . 
At higher densities, dip max /dn decreases and eventually 
becomes negative. The maximum observed rotation was 
ss 10 rad with the applied magnetic field of ~ 0.6 G. 
(This shows the effect of power broadening on the mag- 
netic field dependence; at low light power, the maximum 
rotation would occur at a magnetic field about an order of 
magnitude smaller.) Rotation slope dip/dB\B=o, an im- 
portant parameter for magnetometry, was also measured 
as a function of atomic density in Ref. § . 

The goal of the present contribution is to analyze the 
scaling of optimal magnetometric sensitivity with respect 
to optical thickness of the sample. We analyze NMOR 
for a simple system — an isolated J = 1 — > J = transi- 
tion — for which analytical solutions for the density ma- 
trix are readily obtained. We assume that the transverse 
thickness of the sample is small, so that the trapping of 
spontaneously-emitted radiation can be neglected. (This 
assumption may not be justified in actual experiments 
with optically thick media j9|.) We also neglect mix- 



ing between different velocity groups due to velocity- 
changing collisions. Such mixing is important in the pres- 
ence of a buffer gas, or for anti-relaxation-coated cells. 

A theoretical treatment of NMOR in dense media has 
been carried out earlier in Refs. @ and jTJ. The present 
analysis partially overlaps with this earlier work and ex- 
tends it in several ways: we extend the treatment to 
the case of large Doppler broadening and the interme- 
diate regime where power and Doppler widths are com- 
parable; provide a qualitative discussion of the contribu- 
tion to optical rotation as a function of optical depth in 
the medium; discuss the role of the coherence and the 
alignment-to-orientation conversion effects |Tj|; discuss 
scaling of magnetometric sensitivity with optical thick- 
ness of the medium; and finally, give a general argument 
on scaling based on the observation that an optimized 
nonlinear optical rotation measurement is a way to carry 
out a spin measurement of an atomic ensemble with sensi- 
tivity given by fundamental quantum noise limits |l2| . 

II. DESCRIPTION OF THE DENSITY MATRIX 
CALCULATION 

The calculation is based on a standard density matrix 
approach. The time evolution of the atomic density ma- 
trix p is given by the Liouville equation (see, e.g., Ref. 

where the square brackets denote the commutator and 
the curly brackets the anti-commutator, and the to- 
tal Hamiltonian H is the sum of light-atom interaction 
Hamiltonian Hl = —d ■ E (where E is the electric field 
vector, and d is the electric dipole operator), the mag- 
netic field-atom interaction Hamiltonian Hb = —ft ■ B 
(where B is the magnetic field and pi is the magnetic mo- 
ment), and the unperturbed Hamiltonian Hq. T is the 
relaxation matrix (diagonal in the collision-free approxi- 
mation) 

(ZJm\r\tJm) =7 + 7o<5(£,£ e ), (2) 
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where 7 is the ground state depolarization rate (e.g., due 
to transit of atoms through the laser beam), 70 is the 



2 



spontaneous decay rate from the upper state, and £ rep- 
resents the quantum number distinguishing the ground 
state (£ g ) from the excited state (f e ). A = A + A re r°P is 
the pumping term, where the diagonal matrix 

(£ 9 J s m|A%J ff m> = . ™ - (3) 



describes incoherent ground state pumping (po is the 
atomic density), and 



(£,gJ g m\A repop \£ >g J g m') = 70 ( J 9 , m, 1, q\ J e , m e )( J g , m' , 1, q\ J e , m' e )p^ Jcm ^ c j cm / , (4) 



describes rcpopulation due to spontaneous relaxation 
from the upper level (see, e.g., Ref. Here (. . . | . . . ) 

are the Clebsch-Gordan coefficients. 

The electric field vector is written (see, e.g., Ref. |l6| ) 
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where w is the light frequency, k — lj/c is the vacuum 
wave number, Eq is the electric field amplitude, tp is the 
polarization angle, e is the ellipticity (arctangent of the 
ratio of the major and minor axes of the polarization 
ellipse), and <f> is the overall phase. By substituting (||) 
into the wave equation 
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where P = Tr(pd) is the polarization of the medium, the 
absorption, rotation, phase shift, and change of elliptic- 
ity per unit distance for an optically thin medium can 
be found in terms of the density matrix elements (these 
expressions are given in Ref. fl7[| ). Once the solutions 
for the density matrix are obtained, we will perform an 
integration to generalize the result to media of arbitrary 
thickness. 



III. THE DOPPLER-FREE CASE 

We consider the case ofaJ = l^J = transi- 
tion, and linearly-polarized incident light, with a mag- 
netic field directed along the light propagation direction 
(Faraday geometry). Using the rotating wave approxima- 
tion, the solution of Eq. (Q) is obtained, and from this, 
analytic expressions for thin medium absorption and ro- 
tation are found. These expressions can be simplified by 
assuming that 7 <C 70 ■ We first consider the case where 
the power-broadened line width is much greater than the 



Dopplcr width. In this case, the absorption coefficient 
per unit length a is found to be 



a 



(7) 



(2A/ 7o ) 2 + 2k/3 + I ' 
where A is the light-frequency detuning from resonance 
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is the optical pumping saturation parameter, and 
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is the unsaturated absorption coefficient on resonance, 
where A is the transition wavelength and n is the atomic 
density. For Q <C 7, where Q = gfiB is the Larmor 
frequency (g is the Lande factor, and \x is the Bohr 
magneton) , the slope of optical rotation per unit length, 
dip/(dttdx), (proportional to rotation for small magnetic 
fields) is found to be 
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dndx ~ 7 (2A/ 7o ) 2 + 2k/3' 



(10) 



where we have neglected linear optical rotation. In 
general, optical rotation can be induced by either lin- 
ear dichroism or circular birefringence. Analysis of the 
steady-state polarization of the ground state shows that 
both processes contribute here: the contribution due to 
linear dichroism induced in the medium is given by 



d,p 



dVtdx 
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and the contribution due to circular birefringence (arising 
due to alignment-to-orientation conversion in the pres- 
ence of both the magnetic field and the strong electric 
field of the light |13 ) is given by 
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The sum of the two contributions produces the 
Lorentzian line shape of Eq. (|l^) . 

We now generalize the formulas for absorption and 
rotation to the case of thick media, and find the mag- 
netometric sensitivity. In the Doppler-free case we are 
presently considering, we can further simplify the ex- 
pressions by assuming A = (we will have to include 
non-zero detunings in the discussion of the Doppler- 
broadened case below), and n ^> 1. (We will see from 
the final result that this holds everywhere in the medium 
when the input light power is optimal.) Generally, the 
medium in the presence of a magnetic field produces light 
cllipticity as well as rotation; however, the ellipticity is 
an odd function of detuning and is zero on resonance. 
The rate of change of the saturation parameter as light 
travels through the medium is given by 



dn(x) 
dx 



-0>k(x) 



(13) 



so solving for the saturation parameter as a function of 
position, 



k(x) «Ko - -a x, 



(14) 



where kq is the saturation parameter at x = 0. 

The contribution to the small-field optical rotation of 
the "slice" of the medium at position x is found by sub- 
stituting Eq. (p|) into Eq. Q: 
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This contribution is plotted as a function of position in 
Fig. [j] The plot illustrates that the part of the medium 
near its end contributes significantly more to the overall 
rotation than the part near its beginning. This is because 
light power, and correspondingly, the light broadening of 
the resonance, is lower at the end of the medium [the 
ratio of the rotation per unit length is approximately 
k(£)/kq w 1 — (3/2)a>ol/Ko\. Integrating over the length 
I of the medium gives the total slope: 
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The slope ( |16[ ) is plotted as a function of cto£ in Fig. |2j. 

The photon shot-noise- limited magnctomctric sensitiv- 
ity SB p3] is given in terms of the number N 1 of trans- 
mitted photons per unit time, the slope of rotation with 
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FIG. 1: Normalized contribution to small-field slope (in radi- 
ans) of a slice of the medium of width dx at position x, given 
by Eq. [[5J We have set ko = 1.8a<o^, the value which is later 
seen to produce the greatest magnetometric sensitivity. 




FIG. 2: Normalized slope (in radians) of NMOR as a func- 
tion of optical thickness of the medium [q^ 1 = (3/2)ao^/«o]- 
The value of g _1 that produces the greatest magnetometric 
sensitivity is shown with a dashed line. 



respect to B, and the measurement time t: 



dip 




where A is the cross-sectional area of the light beam and 
we have made the change of variables kq = 3qao(/2. The 
factor 



q 



lln 



1 



reaches a maximum of ~ 0.8 at q « 1.2 so we see that for 
media of sufficient thickness, i.e. where ctol 3> 3, and for 
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the optimum initial saturation parameter kq ~ l.Sctoi, 
we have from Eq. ( |l4[ ) that k(x) ^> 1 for all x. For the 
optimum light intensity, then, 
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This result is consistent with the general observation g{ 
that optimized NMOR provides a method for measuring 
a spin-system at the standard quantum limit (SQL) given 
by spin-projection noise p3| ]. The sensitivity is then ex- 
pected to scale as the square root of the product of the 
number of available atoms and the spin-relaxation time 
(see, e.g., Ref. Sec. 3.1.1), which is indeed the result 



IV. THE DOPPLER-BROADENED CASE 

Now we consider the case where the Doppler width T 
is much greater than the power-broadened line width. In 
this case, we need to average over the over the atomic ve- 
locity distribution, which is equivalent to averaging over 
the Doppler-free spectral profiles. On resonance with the 
Doppler broadened transition, the absorption coefficient 
is given by 
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In this case the Doppler-broadened unsaturated absorb- 
tion coefficient is given in terms of the Doppler-free un- 
saturated absorbtion coefficient by 
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Comparing Eqs. (19) and ([?]), we see that we have re- 
produced a well-known result (see, e.g., Ref. Q, Sec. 
7.2.1) that resonant absorption falls as 1/k for Doppler- 
free media and as l/v^ f° r Doppler-broadened media 
when k> 1. The change in k per unit length is 
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and solving for Ki clS <i> function of position, 
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Note that the behavior of the saturation parameter as a 
function of distance is different in the Doppler-broadened 
case compared to the Doppler-free case [cf. Eq. (|l4|)] 
where the saturation parameter falls approximately lin- 
early with distance. 

Taking the average of small-field rotation over the 
Doppler distribution gives 
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We see that the rotation per unit length scales as 1/ y/n, 
in contrast to the 1/k scaling for the Doppler-free case 
[similar to the situation with absorption, Eqs. (M,^)]. 
This is because the number of atoms producing the effect 
is, in a sense, not fixed; with increasing light power, a 
larger fraction of the Doppler distribution is involved. 

Substituting ( p^ ) into Eq. ( p3| ) to find the contribution 
to the slope as a function of position gives 
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It is interesting to note that while the light power and 
the rotation slope per unit length behave differently in 
the Doppler-broadened case compared to the Doppler- 
free case, Eq. (^4|) has the same functional form as for 
the Doppler-free case [Eq. (|l5|)l. 

Integrating over the length of the medium, we obtain 



dp 

dil DB 
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The behavior of ( Ba ) as a function of a.§l is qualitatively 
similar to that of Eq. ( |l6| ) shown in Fig. |^. However, the 
dependence on kq is different. 

The magnetometric sensitivity is given by 
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with the change of variables 
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The factor 



(p-l)ln 



p-l 



goes to unity as p goes to infinity; it is approximately 0.9 
at p = 5. (Further gain from increased power is minimal, 
and if the power becomes too high, the approximation 
of the Doppler width being much larger than the power- 
broadened width breaks down.) Thus, in this case, 
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for sufficiently high «q. 

This result, where sensitivity increases linearly with 
optical thickness, holds for the case where the power- 
broadened width r~j 70 y k(x) is smaller than the Doppler 
width for all x within the sample, i.e. 



1 « no < (r/70) 2 and 

r/ 7o «a ^«(r/ 7o ) 2 , 



(29) 



since a$l is related to kq by Eq. (|27|). As power and 
optical thickness are increased beyond this range, i.e. 



k > (r/70) and 
» (r/7o) 2 , 



(30) 



we obtain the Doppler-free case, where sensitivity in- 
creases as the square root of the thickness [Eq. (18)]. 



V. THE GENERAL CASE 

A numerical result can be obtained for the general case 
where the restrictions ( p9||3C| ) on k and ao£ are removed. 
In a typical experiment, light power is ~ 1 mW, the laser 
beam diameter is ~ 0.1 cm, A « 800 nm, and r/70 ~ 60. 
Thus the effective ground state relaxation rate due to the 
transit of atoms through the laser beam is 7 « 27T-50 kHz 
and the initial saturation parameter is kq ~ 4 x 10 3 . 
Here, as in a typical experimental procedure, the opti- 
cal depth is varied (by changing atomic density) while 
the laser power is kept constant. Normalized transmis- 
sion, differential small-field rotation, total small-field ro- 
tation, and magnetometric sensitivity are plotted in Fig. 
H as a function of optical depth. For small optical depth, 
k > (r/70) and the medium is effectively Doppler-free. 
Transmission (cx k) falls linearly until the transition to 
the Doppler-broadened case is made (dashed line) . Then 
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FIG. 3: Normalized (a) saturation parameter, (b) differential 
small-field rotation, (c) total small-field rotation, and (d) in- 
verse magnetometric sensitivity as a function of optical depth 
with initial saturation parameter kq — 4 x 10 3 and r/70 = 60. 
Plots (b), (c), and (d) are in units of radians. The dashed line 
indicates the transition from the Doppler-free to the Doppler- 
broadened regime [k = (r/70) 2 ], and the dotted line indicates 
the point at which non-linear effects begin to turn off (« = 1). 
Linear optical rotation is neglected in this plot. The solid line 
indicates the optical depth at which maximum sensitivity is 
achieved. 



transmission falls quadratically until the linear regime is 
reached (dotted line), after which it falls exponentially. 
Differential small-field rotation (oc dip / '[df2d(ao^)] 7 ) ini- 
tially rises, as K falls and power broadening is reduced, 
until non-linear effects begin to turn off. (Linear optical 
rotation is neglected in this plot.) Since magnetomet- 
ric sensitivity depends both on total optical rotation and 
transmission, an intermediate value for the optical depth 
produces the greatest sensitivity (solid line) . Multiplying 
the normalized inverse sensitivity ckp/dQ ■ Jt/k by 



2V2vr 



9V 



h AV7 



(g/ 



(31) 



6 





20 




7 C 
1 J 


1 






10 


1 




3. 


5 









30 








20 






a 


10 










(a) 




_ (ft) 







tion. Maximum rotation is plotted as a function of op- 
tical depth in Fig. [|, for the same parameters as used 
for Fig. ||. The rotation initially rises linearly, the same 
behavior seen in the experiment B mentioned in the in- 
troduction. For large optical depth, however, rotation in 
Fig. |J then begins to rise more quickly, and finally satu- 
rates, whereas experimentally a slower increase and then 
a decrease in rotation is seen. This is evidence for an ad- 
ditional relaxation mechanism not accounted for in the 
present theory; in Ref. || it is attributed to the effect of 
radiation trapping. 
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FIG. 4: (a) Maximum optical rotation and (b) normalized 
magnetic field at which rotation is maximum as a function of 
optical depth. Parameters are the same as those used for Fig. 
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VII. CONCLUSION 



gives the absolute magnitude of sensitivity, ~ 2 x 
1CP 10 G/VHz (we assume g = 1). Although this sen- 
sitivity is not as high as could be achieved with low 
atomic density paraffin-coated cells (~3x 1CP 12 G/VHz 
||), it is, nevertheless, sufficiently high to be of interest 
in practical applications JfO| . In particular, the power- 
broadening of the magnetic field dependence of optical 
rotation at high light power provides an increased dy- 
namic range for magnetometry over the low-power case. 
There are, however, techniques for shifting the narrow 
resonance obtained with a paraffin-coated cell to higher 
magnetic fields, e.g., frequency modulation of the laser 
light @. 

VI. LARGE-FIELD OPTICAL ROTATION 

The maximum optical rotation with respect to mag- 
netic field can be determined with a numerical calcula- 



In conclusion, we have analyzed magnetometric sen- 
sitivity of NMOR measurements optimized with respect 
to light intensity in the case of negligible Doppler broad- 
ening, and in the case of large Doppler broadening. In 
the former case, we find that the sensitivity improves as 
the square root of optical density, while in the latter, it 
improves linearly. In the present discussion, we have ne- 
glected the effect of velocity-changing collisions, which 
makes this analysis not directly applicable to buffer-gas 
and anti-relaxation-coated cells. However, since there is 
full mixing between velocity components in these cells, 
one can expect that the sensitivity should scale as square 
root of optical density (if this quantity can be varied in- 
dependently of the ground state relaxation rate) . 
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